Stabilization mechanism of edge states in graphene 
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It has been known that edge states of a graphite ribbon are zero-energy, localized eigenstates. 
We show that next nearest-neighbor hopping process decreases the energy of the edge states at 
zigzag edge with respect to the Fermi energy. The energy reduction of the edge states is calculated 
analytically by first-order perturbation theory and numerically. The resultant model is consistent 
with the peak of recent scanning tunneling spectroscopy measurements. 



Carbon-based materials have attracted much attention 
from various points of viewi In particular, their electri- 
cal properties are of great interest, where topology plays 
an important role because it is relevant for a rich vari- 
ety of the electronic properties. As one topological as- 
pect, boundaries can induce localized states called edge 
states^* 3 , at graphite edge. Theoretically, edge states are 
zero-energy eigenstates relative to the Fermi energy and 
are predicted to make a certain magnetic ordering^^ Ex- 
perimentally, by scanning tunneling microscopy (STM) 
and spectroscopy (STS) of graphite edge, a peak in the 
local density of states (LDOS) has been observed, and 
it can be identified as the edge statesASi^ An inter- 
esting point is that the peak is located not just at the 
Fermi energy but below the Fermi energy by about 20 
meV. 5 ' 6 Since several possible perturbations would shift 
the energy eigenvalue above the Fermi energy^S it is not 
a simple problem to find a consistent perturbation that 
can decrease or stabilize the energy of the edge state. 
In this letter we show that the next nearest-neighbor 
(NNN) hopping process is a key factor which decreases 
the energy eigenvalue of the edge state. This is shown by 
the first order perturbation theory for the tight-binding 
Hamiltonian and by numerical energy band structure cal- 
culation. 

First, we explain the energy band structures for 
the graphene, using the nearest-neighbor tight-binding 
Hamiltonian for the honeycomb lattice with the hop- 
ping integral —70 (~ —3 eV). We ignore the electron 
spin for simplicity. The eigenenergies of this model are 
given by ±7o|/(k)| where /(k) = £0=1,2,3 e tk-R ° and R a 
(a = 1, 2, 3) are the vectors from an A-sublattice site to 
the neighboring B-sublattice sites. The two energy bands 
are degenerate at the two k points called the K and K' 
points. Now we consider a zigzag nanotube illustrated 
in Fig. ^ Henceforth we define coordinate axes around 
and along the nanotube axis as x\ and x%- In a zigzag 
nanotube (of finite length) the dimcnsionless wave vec- 
tor around the tube, q = k • ai = \/3k\a cc , remains a 
good quantum number and is now quantized as an inte- 
ger multiple of 27r/n with the chiral vector Ch = (n,0). 
For n — ► 00 the model represents a graphite ribbon 
with zigzag edges discussed in Ref. 0- The delocalized 
eigenstates for the nanotube are similar to those of the 



graphene. We take the Brillouin zone to be < q < 2tt. 
The states at the Fermi level without doping are char- 
acterized by /(k) = 0, which implies q = 2n/3, 47r/3, 
corresponding to the K and K' points, respectively. 
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FIG. 1: Lattice structure of a zigzag carbon nanotube with 
finite length. The filled (open) circle indicates the A sublattice 
(B sublattice). Both the left and right edges are zigzag edges, 
ai = \/3fflccei is the primitive vector around tubule axis where 
a cc denotes the carbon-carbon bond length and ei (e2) is the 
unit vector around (along) the tubule axis. 

On the other hand, as we noted in the introduction, 
the zigzag nanotube has edge states in addition to delo- 
calized (bulk) states. Whether or not such edge states are 
allowed depends on the boundary conditions. To analyze 
the edge states, we consider the eigenvalue equation 



-E4> J A 



7o^ +1 

70</>A 



(J = 0,...,iV-l) (1) 



[4>ai 4>b) is t ne wave func- 
the wave function at site 
can be written as 



where G = 270 cos (q/2) and ' 
tion at Jth site (see Fig. [IJ; 
(I, J) of A sublattice (B sublattice, 

^(A,B)l(l) = ex P( iJ ?)<^(A,B)(?)' 

The energy eigenvalues for the nanotube with zigzag 
edges can be calculated by imposing the boundary condi- 
tion: ~E<j)% = G(f>% and —E<j)% = G<P°a- This boundary 
condition for the tube supports edge states around the 
Fermi energy^ By direct calculation one finds that the 
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edge states exist if |G| < 70, implying 2ir/3 < q < Att/3, 
and the wave function of the edge state for —70 < G < 



sinh( J + 1)93 



sinh </3 
70 sinh J<p 
G sinh tp 



sinh( J + l)(f 



sinhtp 



(J = 0,...,N) 



(2) 



where ip is a positive number satisfying 70 sinh(7V + 
1)</? + Gsinh(7V + 2)ip — and depends on the wave 
vector around the tube via the boundary condition: 
e~ v w — 2cos(q/2). The energy eigenvalue is E = 
±70 sinh tp/ sinh(iV + 2)ip, which is exponentially close to 
the Fermi level. When ip is large, the edge state asymp- 
totically behaves exponentially near the edges 



< J< N 
< J < N 



0, 



h N(J-N)<p 



(3) 



where 



=F0^. The localization length is given 
by iy9 _1 (3a cc /2). The energy eigenvalues are E ss 
±7o e_Arv> ~ 0. For < G < 70, the edge state wave 
function is obtained by multiplying (— 1) J by <f> J A and (f>g 
in Eq. © with G — > — G. When approaching G — > ±70, 
i.e., q — > 27r/3,47r/3, the localization length becomes in- 
finite, and the edge state finally becomes a bulk state. In 
this sense, the zero-energy states at q = 27r/3,47r/3 can 
be called critical states^ On the other hand, the state 
with q = 7r corresponds to ip = 00; it is the most local- 
ized state, which has nonzero amplitude only at the edge 
sites. 

Thus far we have considered the nearest-neighbor hop- 
ping process only. Now, we include the NNN hopping 
process with hopping integral — 7„. We show that the 
inclusion of the NNN hopping decreases the edge state 
energy with respect to the Fermi level. 

First, we prove this by means of the first-order pertur- 
bation theory. The Hamiltonian matrix reads as 



H 



-7n(l/(k)| : 



-7o/* (k) 



(4) 



- 7 o/(k) - 7 „(|/(k)| 2 -3) 

The new energy dispersion relation is then given by 

£(k)=± 70 |/(k)|-7„|/(k)| 2 + 3 7 „. (5) 

Henceforth we subtract the constant term 37„ in the en- 
ergy -E'(k). Thus the energy shift due to the NNN hop- 
ping is given by AE = — 7„|/(k)| 2 . The critical states 
have /(k) = 0, which yields AE — 0. The edge states 
appearing in the zigzag nanotubes with long length cor- 
respond to /(k) with complex wave vector along the 
nanotube ki{= k • 62). This naively leads to the result 
that the energy shift, AE, of the edge states due to the 
NNN hopping is zero. Nevertheless, it is not correct, as 
we will see here; the existence of the boundaries is crucial 



for the calculation of energy shift. The resulting shift will 
turn out to be negative, thus stabilizing the edge states. 

Within the first-order perturbation theory, the energy 
shift of the edge state is given by the expectation value 
of the NNN hopping with respect to the edge state. As 
a simple example, we consider the energy shift for the 
most localized edge state (q = 7r). The wave function 
of this edge state is given by tp^j = 4>ai = ( — l) 7 /'\/2n 
and zero otherwise. The energy shift is negative and is 
given by — 27„ cos q — 3j n — — 7«. In the similar way, one 
can evaluate the energy shift for general edge states as a 
function of q. From the asymptotic behavior of the wave 
function given by Eq. © , we evaluate the energy shift to 
the first-order in the NNN hopping as 



A£« 7 „ (-l + G 2 /7o) =7n(2cosg+l). 



(6) 



In particular, for the critical states, it reproduces AE = 
as mentioned previously. Among the localized states, the 
negative energy shift, \AE\ — —AE, is largest for the 
most localized state. 

To confirm these results, we numerically diagonalize 
the tight-binding Hamiltonian for a graphite ribbon with 
zigzag edges. We consider the energy band structure of 
a graphite ribbon with N = 20. In Fig. Ufa), we show 
the energy band structure without the NNN hopping: 
7„ = 0. The edge states form a flatband at E = C 2 
which makes a peak in the LDOS shown in Fig. [2Jc). In 
Fig. IHb), we plot the energy band structure including 
the NNN hopping. We set j n = 0.17^ which shifts 
the Fermi level to Ef = 0.37o>ii The critical states are 
located at Ep{= 37„), and the edge states have lower 
energies by — 7„ < AE < 0, in agreement with Eq. JHJ. 
Therefore, we conclude that the energy minimum at q = 
7r shown in Fig. |2Ib) gives a sharp peak in the LDOS. 
In Fig.[2d) we plot the LDOS at several points near the 
zigzag edge as a function of energy measured from the 
shifted Fermi level. The abovementioned peak is clearly 
seen and this is responsible for the peak in the LDOS 
observed by recent experiments^^ 

It is noted that our model does not include the overlap- 
ping integral (s parameter^) which increases (decreases) 
the conduction (valence) band width. To examine the ef- 
fect of s parameter on the edge states, we performed the 
energy band structure calculation in an extended tight- 
binding framework 12 i and found the similar behavior of 
the energy band structure depicted in Fig. [21 (b). The 
most stable edge state is q — it and the energy eigen- 
value is located below the Fermi level by about -f n . The 
s-parameter does not affect the energy spectrum near the 
Fermi level since the effect of the overlapping integral is 
proportional to the energy of the corresponding states 
measured from the Fermi level. 

Here, we discuss the relationship between our work and 
experimental results?'^'- on a peak of LDOS at graphite 
edges. Niimi et al£ observed a clear peak at a zigzag 
edge and found no such signal at an armchair edge. The 
peak in the LDOS is located below the Fermi energy (de- 
fined by zero bias voltage) by about 20 meV. The inten- 
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(a) 7n = (b) ln = O.I70 
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FIG. 2: Energy band structure of zigzag-edge, nano-graphite 
(a) without the NNN hopping and (b) with the NNN hopping 
(7n = O.I70). Local density of states at some points from the 
edges (0.5 - 2.5 nm) (c) without the NNN hopping and (d) 
with the NNN hopping. In (d), E — is taken as the Fermi 
energy. 

sity of the peak depends on the distance from the edge, 
which can be attributed to the localized nature of the 
edge states. Kobayashi et al& observed a similar be- 
havior. In addition, they observed a peak not only at 
a zigzag edge but also at defect points of an armchair 
edge, while they found no peak at a homogeneous arm- 
chair edge. The important point from our viewpoint of 
this letter is that the above two experimental groups ob- 
served the peak located below the Fermi level by about 
10~ 2 7o- However, this property does not seem to be a 
common property for edge states. By using the STM 
and STS, Klusek et alM found peaks of LDOS in the 
energy range of 20-250 meV above the Fermi level at the 
edges of circular pits on graphite surface. Although the 
circular pits have a mixture of zigzag and armchair edge 
shapes, it is expected that there appear the edge states 
near (local) zigzag edges in the circular pits since Nakada 
et alH showed numerically that localized states appear 
not only in the zigzag edges but also in edges with other 
shapes. Such a general edge state is beyond the scope of 



this letter. It is important to note that samples examined 
by each experimental group were not prepared under the 
same condition. Kobayashi et al& observed the graphcne 
with clear edge structures which were terminated with 
hydrogen in ultrahigh vacuum (UHV) conditions. Their 
procedure of sample preparation can exclude functional 
groups including oxygen which tend to lower the energy 
of the edge states. On the other hand, there is a possi- 
bility that samples used by Niimi et aim* and Klusek et 
al±— included such functional groups since their samples 
were not treated with hydrogen which is activated at high 
temperatures in UHV conditions. 

Although we have demonstrated that the NNN hop- 
ping can decrease the energy of the edge state at zigzag 
edge, there is still a gap between our result, AE — —j n as 
— 10 _1 7o, and experimental data^*& — 10 _2 7o- This gap 
can be attributed to several physical origins. Among 
several factors, the Coulomb interaction would give a 
significant charging energy to the localized edge states. 
The energy shift due to the NNN hopping depends on 
the localization length of the edge state and varies from 
— In < AE < 0. The most localized edge states have the 
largest energy shift AE as — 7„; as the localization length 
becomes longer, the energy shift approaches AE as 0. Be- 
cause the Coulomb charging energy is basically inversely 
proportional to the localization length, the energy shift 
of the edge state relative to the Fermi level will be re- 
duced. It is noted that the Coulomb interaction is of 
particular importance from the point of view of the spin 
polarization^ 

In summary, we point out that the NNN hopping pro- 
cess decreases the energy of edge states at zigzag edges 
with respect to the Fermi energy. The energy reduction 
depends on the localization length of the edge states. 
The most localized edge states have the largest energy 
reduction due to the NNN hopping, while the critical 
states stay on the Fermi level when the sample is suf- 
ficiently large. We calculate the energy shift [Eq. 
by the first-order perturbation theory, and confirm the 
result numerically as shown in Fig. [21 
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